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We make a theoretical study of the quantum oscillations of the tunneling magnetoresistance 
(TMR) as a function of the spacer layer thickness. Such oscillations were recently observed in 
tunneling junctions with a nonmagnetic metallic spacer at the barrier-electrode interface. It is 
shown that momentum selection due to the insulating barrier and conduction via quantum well 
states in the spacer, mediated by diffusive scattering caused by disorder, are essential features 
required to explain the observed period of oscillation in the TMR ratio and its asymptotic value for 
thick nonmagnetic spacer. 

PACS numbers: 75.70.-i, 75.70.Pa, 73.40.Gk 



Large magnetoresistance jl], || observed in ferromag- 
netic tunneling junctions such as Fe/Al 2 03/Fe and 
Co/A^Oa/CoFe currently attracts much interest due to 
the possibility of its application to magnetic sensors and 
MRAM elements. Because the tunneling magnetoresis- 
tance (TMR) ratio is related to the spin polarization of 
the ferromagnetic leads || ||, attempts have been made 
to fabricate junctions with more highly spin-polarized 
ferromagnets ]|, ||. Realistic calculations @, ||, on the 
other hand, have given much higher TMR ratios than the 
observed values, which is probably due to their assump- 
tion of epitaxial structures. Recent experiments on TMR 
using epitaxial junctions 10, [H], however, were un- 



successful in producing TMR ratios as high as expected. 
Thus, our understanding of the relationship between the 
electronic structure of the ferromagnets and the TMR 
ratio is far from complete. 

The most important factor governing the TMR ra- 
tio may be the electronic structure at junction inter- 
faces 
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In order to clarify its role, several ex- 
periments have been performed to measure the depen- 
dence of TMR ratio on the thickness of a nonmagnetic 
metal layer inserted at the interface JlJ, EM [l(|. The 
observed TMR ratios show almost monotonic decrease 
with increasing thicknesses of inserted layers of Au, Cu, 
or Cr, contrary to a theoretical study for clean junctions 
jl7j which shows clear oscillations of the TMR ratio as 
a function of the nonmagnetic layer thickness. Zhang 
and Levy Q have successfully explained this decrease 
in TMR ratio in terms of the decoherence of electron 
propagation across a nonmagnetic layer. However, re- 
cent experiments by Yuasa et al. show clear oscillations 
of the TMR ratio as a function of Cu layer thickness 
for high quality NiFe/A^Os/Cu/Co junctions in which 
the Co/Cu electrode is a single crystal In their ex- 
periments, two characteristic features of the oscillations 



have been observed: (i) the average TMR ratio decays 
to zero with increasing nonmagnetic layer thickness; (ii) 
the period of the oscillations is determined solely by the 
belly or long period Fermi wave vector fcp of Cu. The 
observed period agrees quite well with that of the oscil- 
lations of photoemission spectra caused by quantum well 
states in Co/Cu multilayers [^(J. From the theoretical 
point of view this is confusing, since in addition to the 
Fermi wave vector |2l|] , another wave vector, i.e., the cut- 
off k-point k cp , given by the depth of the quantum well, is 
also known to contribute to the conductance oscillations 
p2| . This wave vector dominates the predicted oscilla- 
tions of CPP-GMR in a Co/Cu/Co trilayer ||. In fact, 
the calculated oscillations of TMR for a clean junction 
[jlTj cannot be explained by a single period determined 
by only. Furthermore, the asymptotic value of the 
TMR ratio calculated for a thick spacer layer is finite, 
which disagrees with the observed results. The purpose 
of the present work is to reconcile the theoretical results 
with the observed ones and thus deepen our understand- 
ing of the TMR effect. 

In this Letter we will show that the combined effects of 
barrier thickness and disorder can explain the experimen- 
tal results. In particular, we will demonstrate that i) in- 
creasing barrier thickness increases the amplitude of the 
kp oscillation period relative to the k cp oscillation period, 
ii) disorder introduced in the barrier also weakens the am- 
plitude of the fc cp oscillation period, and iii) the disorder 
decreases the asymptotic value of the TMR ratio. These 
results are interpreted in terms of the momentum selec- 
tion of electrons incident on the barrier interface and in 
terms of the diffusive scattering due to disorder which 
opens additional conduction channels via quantum well 
states. In the first part of this Letter, these effects will 
be demonstrated by numerical calculation for a single- 
orbital tight-binding model. We will then demonstrate 
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FIG. 1: Potential profile of a FM/I/NM/FM junction. 



that the calculated results can be reproduced by the sta- 
tionary phase approximation |2j| . This implies that this 
technique is applicable to a realistic multi-orbital tunnel- 
ing junction, where a purely numerical calculation would 
be unfeasible. The results indicate that k? of Cu is really 
responsible for the oscillation period observed. 

Let us consider a FM/1/NM/FM junction on a simple 
cubic lattice with lattice spacing a, where FM, I, and NM 
denote a ferromagnetic electrode, an insulating barrier, 
and a nonmagnetic metallic spacer, respectively. Initially 
we adopt a single-orbital tight-binding Hamiltonian in 
order to model a C0/AI2O3/CU/C0 junction: 



(1) 



where c icr (cj cr ) is the annihilation (creation) operator of 
an electron with spin a at site i, t the hopping integral 
between nearest neighbor sites, and Vi„ the on-site po- 
tential for an electron with spin a at site i. Since the 
majority (+) spin band of Co is similar to the Cu band, 
we assume that Vfm+ is equal to Vnm- Figure 1 shows 
the potential profile of the system. Quantum well states 
are formed in NM only for electrons with minority (— ) 
spin in the right FM. Since the insulating AI2O3 barrier 
is amorphous in real junctions, we introduce disorder in 
the barrier by requiring that Vi a takes Vi+AV or V\ — AV 
values randomly depending on the site in the barrier. 

The Kubo formula and a recursive Green's function 
method are used to calculate the tunneling conductances 

G ++ , G__, G^ , and G |_, where G++ and G are 

the conductances in parallel alignment for \ and J.-spin 

electrons, respectively, and G^ and G \- are those in 

antiparallel alignment for f and j-spin electrons, respec- 
tively. The conductance is given by 
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where £ CTCr '(k|i — > k'||) is the transmission coefficient for 
an electron incident from the left FM with k|| and scat- 
tered to the right FM with k'||. TMR ratio is evaluated 
from the conductances in the parallel and antiparallel 
alignments as TMR = 1 - (G+_ + G_+)/(G ++ + G__ ). 
In order to treat the disorder introduced in the insulating 



barrier, we use the single-site coherent potential approxi- 
mation (CPA). The vertex correction to the conductance, 
which describes diffusive scattering, is calculated consis- 
tently with the coherent potential (self-energy) so that 
the current conservation is satisfied pi]. We have also 
performed numerical simulations for finite-size clus- 
ters and checked that the results obtained by the two 
methods agree. 

In these numerical calculations, we use Vfm+ = Vnm = 
2.382£, V ¥ m- = 5.382£, and Fermi energy E F = 0.0. The 
choice of these parameters gives commensurate periods 
of oscillation as shown below. As for the insulating bar- 
rier, parameters V\ = 9.0t and AV = are used for clean 
junctions, and V\ = 9.0i and AV = 0.5t are used for dis- 
ordered junctions. We only show the calculated results 
for disorder within the insulating barrier. However, we 
have checked that the results are not changed qualita- 
tively even when we introduce disorder at the interface 
between the nonmagnetic spacer and the ferromagnetic 
electrode. 

Figures 2(a) and 2(b) show, respectively, the spin- 
dependent conductances and TMR ratios of junctions 

without disorder. It can be seen that G__ and G^ 

oscillate with the NM layer thickness Lnm due to inter- 
ference effects caused by the quantum well. These os- 
cillations show more than one period. The Fermi wave 
vector kp and the cut-off k-point k cp of NM are given 
by 2icos(fcpa) = Vnm — Ep — 4t and 2tcos(k cp a) = 
Vfm- — Vnm, respectively p2| . Therefore, the periods of 
oscillation estimated from k-p = 47r/5a and fc cp = 2-7r/3a 
are 5a and 3a, respectively. The periods of oscillation 
in G and G^ shown in Fig. 2(a) may be inter- 
preted as a superposition of these two periods as dis- 
cussed later. The situation is analogous to that of CPP- 
GMR in a Co/Cu/Co trilayer Q . The TMR ratio shown 
in Fig. 2(b) oscillates with the same periods as the con- 
ductance and has a finite asymptotic value for large NM 
thicknesses. These results are consistent with the previ- 
ous results M where the effect of disorder was ignored. 

Figure 3 shows the dependence of oscillations in G 

on the barrier thickness L\. Here the conductances are 
normalized to the asymptotic values Gf°_ obtained for 
Lnm — > 00. It can be seen that the oscillation period 
tends to 5a with increasing barrier thickness. This re- 
sult is explained as follows. The wave vector k|| parallel 
to the interface is conserved in the system without dis- 
order, that is, W(k|| -> k'||) = t tro -/(k||)o*k || ,k'|| ■ The 
transmission coefficient depends strongly on the angle of 
incidence of electrons tunneling across the barrier, and 
the normal incidence contributes most to the conduc- 
tance. It follows that, as the barrier thickness increases, 
the oscillation given by cut-off k-points, i.e., k|| 7^ 0, be- 
comes progressively weakened compared to that given by 

the Fermi wave vector of NM, i.e., k|| = 0. As for G^ , 

we could not see the increase in oscillation period from 3a 
to 5a unless we increase L\ further. This might be due 
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FIG. 3: Conductance G__ for various thicknesses of clean in- 
sulating barrier. Conductances are normalized by the asymp- 
totic value G?°_ obtained within the limit of large spacer 
thickness. Dashed lines are visual aids. 



to the fact that L| (k|| = 0)/^ k t_| (k||) is smaller 

than i__(k|| = 0)/^ k i__(k||): that is, the contribu- 
tion of normal incidence to the conductance in G_| is 

less than that in G As a result, the oscillation period 

of the TMR ratio is not quite 5a for the present barrier 
thickness. 

We now introduce disorder into the insulating barrier 
and show the calculated results of the spin-dependent 
conductances and TMR ratios in Figs. 4(a) and 4(b), 
respectively. It can be seen that the conductance G + _ 
is enhanced by disorder whereas the other conductances 
G++, G__, and G |_ are hardly affected. G_| now oscil- 
lates almost exclusively with period 5a, (i.e. fcp period) 
about G ++ (see Fig. 4(a)). This results in a TMR ra- 
tio which is decreased and oscillates around zero with 
period 5a. This should be contrasted with the ordered 



FIG. 4: Spin dependent conductance (a) and TMR ratio (b) 
calculated for disordered junctions. Conductances for j and 
I spin electrons are plotted by triangles and circles, respec- 
tively, for parallel (solid symbols) and antiparallel alignment 
(open symbols) of magnetizations. Solid lines indicate results 
obtained by the stationary phase approximation while dashed 
lines are visual aids. Inset of (b): Asymptotic values of TMR 
ratio TMR°° obtained within the limit of large spacer thick- 
ness calculated for clean (open circles) and disordered (solid 
circles) junctions. 



case in which the TMR ratio oscillates with a mixed pe- 
riod about a constant background (cf. Fig. 4(b) and 
Fig. 2(b)). The asymptotic values of the TMR ratio as 
^nm - * oo are shown in the inset of Fig. 4(b) as func- 
tions of the barrier thickness. Both the cases with and 
without disorder are shown. The asymptotic value of 
the TMR ratio of junctions without disorder decreases 
slowly with increasing Lj, whereas that of junctions with 
disorder decreases rapidly and becomes zero for large L\. 

To gain a better understanding of the effects of dis- 
order on the magnitude of G^ , and on the period 

of oscillations, we have calculated the dependence of 
the transmission coefficient on k||. Figures 5(a) and 
5(b) show the transmission coefficients (kii), where 



T^ky) = X) k . t 



k'li) is the transmission co- 



efficient for an electron incident from the left FM with 
momentum kn on the barrier without and with disorder, 
respectively. When there is no disorder, the contribu- 
tion to T-\ in the momentum space is concentrated near 

k|| = (0,0). However, inclusion of disorder gives rise to 

additional contributions to of momenta outside this 

area. This is due to the fact that k|| need not be con- 
served in diffusive scattering. In the absence of disorder, 
only ky points on the Fermi surfaces, that satisfy the 
k|| conservation, may contribute to the conductance. For 
diffusive scattering, on the other hand, the entire set of ky 
points on the Fermi surface contributes to the transport. 
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FIG. 5: ky-dependence of transmission coefficients T+_ cal- 
culated for (a) clean and (b) disordered junctions. 



More precisely, k|| points corresponding to these quan- 
tum well states contribute to the conductance. These ky 
points appear as spiky peaks in Fig. 5(b) and they fall 
on concentric rings in the kii space. 

It is clear from Fig. 5(b) that the number of open kii 
channels contributing to T-\ , is the same as that con- 
tributing to T ++ . This explains the increase in the con- 
stant part of the conductance G_| to a value of approx- 
imately G++- In addition, the introduction of diffusive 
scattering has almost eliminated the sharp momentum 
cut-off observed in fig. 5(a), which explains why the fc cp 

oscillation period of G^ is weakened by disorder. The 

other transmission coefficients are not greatly affected by 
the introduction of disorder, as scattering cannot open 
new k|| channels for these cases. This explains why the 
introduction of disorder has little affect on the conduc- 
tances G++, G__, and G |_. 

We therefore expect that in the presence of disorder, 
the oscillatory part of the conductance is derived entirely 
from states in the region of ky = (0, 0). In order to check 
this hypothesis, we use the stationary phase method (see 
Ref. j2^| and references therein), which is able to deter- 
mine the oscillatory contributions from isolated regions 
of the Brillouin zone, for thick spacers. The results, de- 
picted by solid curves in fig. 4(a) and 4(b), are in excellent 



agreement with the numerical calculations for Lnm 5. 
This fact indicates that the oscillation period observed 
in the experiments may be determined in the stationary 
phase approximation for realistic systems. The experi- 
mental finding jl9| that the oscillation period is deter- 
mined by &f of Cu spacer is thus naturally explained. 
Realistic calculation for the TMR oscillation and its bias 
dependence is in progress. 

In summary, the period of oscillation determined by k-p 
of the spacer is dominant in TMR due to the ky -selection 
by the insulating barrier and the deconfinement of the 
quantum well states by disorder. The diffusive scatter- 
ing caused by the disorder increases the conductance in 
antiparallel alignment by opening new conductance chan- 
nels via quantum well states and results in the oscillation 
of the TMR ratio around an averaged value close to zero. 
The success of the stationary phase approximation in re- 
producing the numerical results indicates that the oscil- 
lation period observed in realistic tunneling junctions can 
be explained in terms of fcp of Cu spacer. 
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